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I.  INTRODUCTION 

Natural  draft  hyperbolic  cooling  towers  have  been  used, 
in  this  country,  for  over  two  decades.  They  have  increased 
in  height  during  this  time  by  over  twenty  percent  to  the 
present  height  of  550  feet  (170  meters).  Towers  of  650  feet 
(200  meters)  and  more,  in  height,  are  being  considered  for 
future  construction.  While  the  typical  shell  thickness  is 
approximately  15  inches  (0.38  meters),  the  typical  base 
diameter  is  300  feet  (92.3  meters). 

While  their  unique  shape  is  both  structurally  and 
thermally  efficient,  they  are  subject  to  dynamic  loading, 
especially  from  wind. 

A  considerable  amount  of  research  has  been  conducted  on 
hyperbolic  cooling  towers.  This  research  includes  the  use  of 
analytical  (mathematical),  small-scale,  and  full-scale  pro- 
totypical models. 

An  analytical  model  consists  of  mathematical  formula- 
tions that  are  used  to  calculate  or  describe  various  proper- 
ties and  responses  of  a  structural  system.  There  are  many 
problems  associated  with  the  utilization  of  a  purely  theo- 
retical model.  The  accuracy  of  such  a  model  is  dependent 
upon  various  assumptions  made  by  the  user.  These  may  in- 
clude stiffness  properties,  boundary  conditions,  and  sus- 
pected foundation/soil  interactions. 

There  are  also  problems  associated  with  the  utilization 
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of  small-scale  models.  It  is  very  difficult  to  construct  a 
model  with  all  of  the  required  parameters.  These  may  include 
Euler  numbers,  as  well  as  mass  density,  Young's  modulus, 
Poisson's  ratio,  and  roughness  coefficients. 

Scale  models  of  reinforced  concrete  present  some  spec- 
ial problems.  For  example,  full-scale  reinforcing  steel  is 
deformed,  or  roughened,  to  enhance  the  bond  between  the 
concrete  and  steel.  The  effect  of  this  enhancement  is  very 
difficult  to  reproduce. 

In  the  case  of  a  hyperbolic  cooling  tower,  the  geometry 
of  the  structure  may  also  present  some  special  difficulties. 
The  shell  thickness  of  a  scale  model  may  preclude  the  use  of 
concrete  as  the  modeling  medium. 

The  testing  of  a  full-scale  structure  presents  some 
unique  problems.  It  is  not  economically  feasible  to  conduct 
destructive  tests  on  a  full-scale  structure.  Most  non- 
destructive tests  are  employed  to  determine  in-situ  material 
properties  (swiss  hammer,  ultra-sound),  or  connection  integ- 
rity (x-ray).  The  physical  size  of  the  structure  may  present 
problems  with  accessibility  to  the  researcher,  and  often, 
information  gained  from  the  testing  of  one  structure  may  not 
be  applicable  to  other  structures  of  the  same  type  (13). 

In  general,  the  experimental  research  done  on  hyper- 
bolic cooling  towers  has  focused  on  defining  the  differen- 
tial wind  and  pressure  loadings.  One  exception  is  the 


research  by  Abel  and  Billington,  who  compared  the  dynamic 
behavior  of  cracked  and  uncracked  hyperbolic  shells  (4) . 
Another  is  the  development  by  Lee  and  Gould  of  a  complex 
algorithm  for  the  dynamic  analysis  of  thin  shell  structures 
supported  on  interactive  foundations  (15). 

Other  research  includes  a  study  to  collect  and  inter- 
pret the  results  of  full-scale  wind  pressure  measurements, 
and  to  determine  the  corresponding  displacements  and  forces 
through  finite  element  analysis,  by  Basu  and  Gould  (7). 

Sollenberger ,  et.al.,  established  the  level  and  dis- 
tribution of  design  dynamic  wind  pressure  loadings  over  a 
typical  cooling  tower  (19),  and  supported  their  conclusions 
with  data  obtained  from  experimental  measurements  on  the 
tower  at  Martin's  Creek,  PA.  Random  wind  response  of  a 
cooling  tower  was  performed  by  Kapania  and  Yang,  using  a 
Monte  Carlo  simulation  approach  (14) .  Results  were  presented 
in  the  form  of  a  time  history  response,  and  compared  to 
those  obtained  in  the  frequency  domain. 

Abu-Sitta  and  Hashish  have  presented  an  approach  to 
predict  dynamic  wind  stresses  in  upright  shells  of  revolu- 
tion, and  compared  those  predictions  to  the  measured  re- 
sponse of  an  aero-elastic  wind  tunnel  model  (5). 

Scale  models,  on  the  order  of  1/250  and  1/500,  of 
multiple  tower  groups,  including  surrounding  structures  and 
topography,  were  constructed  and  wind  tunnel  tested  by 


Armitt  (6).  A  1/40  scale  micro-concrete  hyperbolic  tower 
was  constructed  by  Swartz ,  et.al.,  and  subjected  to  an 
internal  vacuum  to  simulate  external  pressure  (21).  The 
tower  was  tested  to  buckling,  repaired,  and  tested  again. 
A  second  tower  is  presently  under  construction. 

Other  studies  include  a  post  mortem  on  the  towers  at 
Ferrybridge,  England  (17),  and  Ardeer,  Scotland  (18),  as 
well  as  a  study  of  the  damage  inflicted  by  a  mobile  crane, 
in  the  wake  of  a  tornado,  and  the  subsequent  repair  thereof, 
on  the  Grand  Gulf  cooling  tower,  in  Mississippi  (10). 

A  better  understanding  of  a  complex  system,  such  as  a 
hyperbolic  cooling  tower,  may  be  enhanced  by  the  utilization 
of  an  analytical  model  in  conjunction  with  the  testing  of  a 
prototype  full-scale  structure. 

One  method  of  full-scale  structural  testing  is  the 
analysis  of  ambient  vibration  data.  Ambient  vibration  data 
analysis  is  used  to  determine  the  natural  frequencies  and 
associated  mode  shapes  of  a  full-scale  structure.  The  term 
'ambient  vibration'  refers  to  culturally  induced  vibrations, 
as  well  as  those  naturally  induced  by  wind,  wave,  and 
micro-seismic  activity. 

Advances  in  the  technology  required  for  ambient  vibra- 
tion data  collection  and  analysis  (accelerometers ,  digital 
computers,  recorders,  etc.)  have  increased  the  feasibility 
of  this  type  of  full-scale  structural  testing.  It  is  rela- 


tively  inexpensive,  does  not  result  in  the  destruction  of 
the  structure,  and  yields  results  based  on  actual  behavior 
of  the  structure  at  low  levels  of  excitation. 

Although  many  types  of  structures  have  been  the  subject 
of  ambient  vibration  studies,  as  yet  (1985),  none  of  these 
tests  have  been  performed  on  hyperbolic  cooling  towers. 
Several  multi-story  buildings  have  been  tested,  including 
the  24-story  steel  frame  office  tower  tested  by  Craig, 
et.al.,  (8).  Thus  structure  was  tested  at  various  stages  of 
construction,  and  the  changes  in  modal  parameters  were  det- 
ermined. Other  tested  multi-story  buildings  include  the 
Imperial  County  Services  Building  (16),  as  well  as  an  18 
story  pre-cast  building  (20). 

Tested  structures  also  include  the  Golden  Gate  Bridge, 
both  the  suspended  structure  (1),  and  the  pier  tower  struc- 
ture (2),  as  well  as  a  400  ft.  long,  symetrical ,  reinforced 
concrete  box  girder  bridge  (9) . 

Other  types  of  structures  subjected  to  ambient  vibra- 
tion analysis  include  dams  (3),  inflatable  structures  (11), 
and  cylindrical  fluid  storage  tanks  (12). 

1.1  Research  Objective 

It  is  the  aim  of  this  research  to  determine  the  dynamic 
characteristics  (natural  frequencies  and  associated  mode 
shapes)  of  a  hyperbolic  cooling  tower. 

Work  on  an  analytical  model  (finite  element  analysis) 
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of  the  tower  is  currently  in  progress,  at  Kansas  State  Univ- 
ersity, and  findings  will  be  reported  to  the  National  Sci- 
ence Foundation. 

It  is  hoped  that  the  information  gained  from  this 
research  will  be  useful  in  calibrating  such  an  analytical 
model,  as  well  as  in  verifying  the  results  of  small-scale 
model  testing  in  the  future. 

For  this  purpose,  a  series  of  ambient  vibration  tests 
was  conducted  on  the  hyperbolic  cooling  tower  at  Portland 
General  Electric 's  Trojan  Nuclear  Generating  Station,  near 
Rainier,  Oregon. 

1.2  Description  of  Structure 

Completed  in  May  of  1973,  the  tower  is  an  offset  shell 
of  revolution  (hyperboloid  of  one  sheet) ,  the  geometry  of 
which  is  given  by  the  following  equation: 
[(R-C)2/a2]-[Z2/b2]=l, 
where  R  is  the  radius  of  revolution  at  height  Z,  and  the 
curve  parameters  a  and  b  equal  34.6  ft.  (10.6  meters),  and 
128  ft.  (39.3  meters)  respectively.  The  offset  distance 
between  the  axis  of  revolution  and  the  centroidal  axis.C, 
is  equal  to  81.4  ft.  (24.8  meters). 

The  tower  also  has  the  following  dimensions: 

Base  Diameter  385.0  ft.  (117.35  m.), 

Lintel  Diameter  364.5  ft.  (111.10  m.), 

Throat  Diameter  232.0  ft.  (  70.72  m.), 
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Cornice  Diameter  249.8  ft.  (  76.14  m.) 

Lintel  Elevation  41.0  ft.  (  12.50  m.). 

Throat  Elevation  393.6  ft.  (119.19  m.). 

Cornice  Elevation  492.0  ft.  (149.97  m.) 

Lintel  Thickness  3.75  ft.  (   1.13  m.), 

and  Throat  Thickness  0.83  ft.  (   0.25  m.). 

A  4.5  ft.  (1.37  en.)  cantilevered  slab  serves  as  a 
stiffening  ring  at  the  cornice  elevation,  and  is  also  used 
as  an  access  walkway. 

Access  ladders  are  located  circumferentially  along  the 
tower,  and  spaced  at  ninety  degree  intervals.  Three  of  these 
ladders  reach  the  420  ft.  (128  m.)  level,  while  the  fourth 
provides  access  to  the  top  of  the  tower. 

Landings  are  located  just  above  the  lintel  beam,  and  at 
approximate  elevations  of  150  ft.  (46tn.),  270  ft.  (82  m.), 
and  390  ft.  (120  m.)  on  all  ladders.  An  additional  landing 
is  provided  on  the  fourth  ladder  at  approximately  485  ft. 
(148  m.).  The  structure  has  an  overall  height  of  499  ft. 
(152  m.). 


II.  EQUIPMENT  AND  PROCEDURE 

Ambient  vibration  data  were  collected  by  placing  acc- 
eleration sensitive  transducers  at  various  locations  on  the 
tower,  and  recording  the  resulting  acceleration  time  history 
on  magnetic  tape  for  later  analysis. 

Accelerometer  locations  were  selected  as  part  of  a 
comprehensive  testing  routine  to  insure  a  sufficient  quant- 
ity of  data  for  the  identification  of  principal  mode  shapes. 

Both  the  equipment  and  test  routine  will  be  discussed 
in  this  chapter. 

2.1  Equipment 

Six  basic  pieces  of  equipment  were  used  to  collect  and 
analyse  the  ambient  vibration  data. 

Accelerometer 

An  accelerometer  is  a  transducer  used  for  the  detection 
and/or  measurement  of  accelerations.  The  transducers  used  in 
this  research  were  FBA-11  Force  Balance  Accelerometers,  man- 
ufactured by  Kinemetrics  Inc.  of  Pasadena,  CA. 

The  FBA-11  is  a  single  axis  instrument  with  a  nominal 
output  of  2.5  volts/g.,  a  nominal  range  of  ±  1.0  g. ,  and  a 
natural  frequency  of  50  Hertz.  They  are  damped  to  70  percent 
of  critical  damping.  Four  such  accelerometers  were  used  in 
the  data  acquisition  system. 

A  modified  C-type  clamp  was  used  to  attach  each  trans- 
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ducer  to  the  tower.  The  clamps  were  modified  to  allow  place- 
ment of  the  accelerometer  in  a  level  or  vertical  position, 
regardless  of  the  position  of  the  clamp  itself.  This  was 
accomplished  by  utilizing  two  pin-type  connections,  the 
axes  of  which  form  right  angles,  between  the  transducer 
and  the  clamp.  This  gave  the  effect  of  a  'universal'  joint. 
Figure  2.1  is  an  exploded  view  of  the  modified  clamp. 

Strong  Motion  Accelerograph 
The  accelerometers  were  powered  by  the  SMA-3  Strong 
Motion  Accelerograph,  also  manufactured  by  Kinemetrics  Inc. 
The  SMA-3  also  provided  various  functional  test  signals. 
These  include  an  operational  test,  as  well  as  a  calibration 
test. 

The  SMA-3  contains  two  12-volt  batteries  and  an  inter- 
nal charging  system.  These  enable  it  to  be  used  in  remote 
testing  situations. 

Signal  Conditioner 
Signal  conditioning  was  provided  by  the  SC-1  Signal 
Conditioner,  manufactured  by  Kinemetrics  Inc.  The  SC-1  is  a 
portable  seismic  instrument  with  four  separate  data  acquis- 
ition channels  providing  adjustable  low-pass  filtering  and 
amplification.  The  SC-1  was  designed  specifically  for  use 
in  data  acquisition  systems  where  adjustable  filtering  or 
amplification  are  desireable. 
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Chart  Recorder 
The  output  from  each  channel  of  conditioned  signal  was 
recorded  on  a  Gould  Brush  Model  4200  four  channel  strip 
chart  recorder.  This  was  done  in  an  effort  to  provide  a 
record  for  on-site  visual  inspection  of  the  response  data. 

Magnetic  Tape  Recorder 
The  output  from  each  channel  was  also  recorded  on  a 
four  channel  Bruel  &  Kjaer  7006  FM  recorder.  This  created 
a  permanent  record  for  future  analysis. 

Spectrum  Analyser 

The  analysis  of  ambient  vibration  signals  was  parti- 
ally accomplished  through  utilization  of  the  Hewlett  Packard 
HP-3582A  real  time,  dual  channel  spectrum  analyser.  The 
analyser  was  used  to  digitize  the  acceleration  signals,  and 
to  compute  the  Fast  Fourier  Transform,  as  well  as  to  deter- 
mine the  coherence  function  and  relative  phase  between  sta- 
tions. These  processes  will  be  explained  further  in  section 
2.3. 

Figure  2.2  is  a  schematic  diagram  of  the  data  acquis- 
ition and  analysis  system. 

2.2  Test  Procedure 

The  collection  of  ambient  vibration  data  was  accomp- 
lished by  placing  the  accelerometers  at  various  locations 
and  orientations  on  the  tower,  and  recording  the  resulting 
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acceleration  time  history  signals  for  later  analysis. 

The  axisymetry  of  the  structure  required  the  formula- 
tion of  a  coordinate  system  for  use  in  locating  the  points 
of  interest  on  the  tower.  The  system  contains  three  para- 
meters. These  are  the  elevation,  position,  and  orientation 
of  the  accelerometer . 

The  base  of  the  tower  was  designated  as  level  0,  and 
each  subsequent  landing  was  assigned  a  level  number. 

The  ladder  which  reached  the  top  of  the  structure  was 
designated  as  position  0,  and  each  subsequent  ladder,  pro- 
ceeding in  a  counter-clockwise  direction,  was  given  a 
position  number  (90,  180,  or  270). 

The  final  parameter  is  the  orientation  of  the  accelero- 
meter. The  accelerometer  was  placed  in  one  of  three  orien- 
tations ; 

1.  Vertical,  pointing  upward  and  designated  V, 

2.  Tangential  to  the  level  circular  plane  of  the 
tower,  pointing  counter-clockwise,  and  desig- 
nated T,  or 

3.  Normal  to  the  axis  of  revolution,  pointing 
inward,  and  designated  N. 

Thus  the  designation  5-0-N  identifies  the  level, 
position,  and  orientation  as  landing  number  5,  ladder  number 
Q ,  and  oriented  normal  to  the  axis  of  revolution,  pointing 
inward. 

Figure  2.3  shows  the  designation  of  parameters. 
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Because  the  amplitude  of  the  tower's  response  is  ex- 
pected to  change  between  tests,  it  is  neccessary  to  normal- 
ize the  response  at  each  location  to  the  response  of  a  single 
continuously  mounted  accelerometer .  Therefore,  one  acceler- 
ometer  was  positioned  for  the  duration  of  the  testing  at  the 
so  called  reference  station.  This  station  provides  a  datum, 
or  point  of  reference, for  the  comparison  of  relative  response 
of  the  other  stations.  The  reference  station  coordinates, 
in  this  case,  were  5-0-N. 

The  order  in  which  the  tests  were  carried  out,  the 
classification  of  each  test,  and  the  location  of  each  trans- 
ducer are  listed  in  Table  2.1. 

This  sequence  of  testing  was  selected  to  contain  the 
desired  test  configurations,  and  to  minimize'  the  number  of 
times  that  each  transducer  had  to  be  relocated.  This  reduced 
the  potential  for  damage  to  the  system,  as  well  as   reducing 
the  effort  required  to  relocate  the  accelerometers. 

The  first  test  was  used  to  calibrate  the  accelerometers. 
This  was  accomplished  by  placing  all  of  the  transducers  at 
one  location  and  recording,  simultaneously,  each  signal  on 
the  chart  recorder.  The  strip  chart  was  then  observed  for 
any  discrepancies.  If  significant  discrepancies  were  noted, 
a  calibration  factor  would  have  been  determined  to  correct 
for  the  discrepancy. 
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Table  2.  1  Test  Routine 


Test   # 

Ch.   1 

Ch.    2 

Ch.   3 

Ch.    4 

Class.   * 

0. 

5-0-N 

5-0-N 

5-0-N 

5-0-N 

Calib. 

1. 

5-0-N 

5-0-N 

5-0-V 

5-0-T 

A 

2. 

5-0-N 

4-0-N 

4-0-V 

4-0-T 

A 

3. 

5-0-N 

4-0-rl 

3%-0-N 

3-0-N 

B 

4. 

5-0-N 

3-0-N 

3-0-V 

3-0-T 

A 

5. 

5-0-N 

3-0-N 

2%-0-N 

2-0-N 

B 

6. 

5-0-N 

2-0-N 

2-0-V 

2-0-T 

A 

r. 

5-0-N 

2-0-N 

1^-0-N 

1-0-N 

B 

8. 

5-0-N 

1-0-N 

1-0-V 

1-0-T 

A 

9. 

5-0-N 

0-0-N 

0-0-V 

0-0-T 

A 

10. 

5-0-N 

2-0-N 

2-90-N 

2-180-N 

C 

11.' 

5-0-N 

2-0-N 

2-90-T 

2-180-N 

C 

12. 

5-0-N 

2-0-V 

2-90-V 

2-180-V 

D 

13. 

5-0-N 

1-0-N 

1-90-N 

1-180-N 

C 

14. 

5-0-N 

1-0-V 

1-90-T 

1-180-N 

A 

15. 

5-0-N 

1-0-V 

1-90-V 

1-180-V 

D 

16. 

5-0-N 

0-0-N 

0-90-N 

0-180-N 

C 

17. 

5-0-N 

0-0-V 

0-90-V 

0-180-V 

D 

*  A=Triaxial  Response 
B=Lateral  Mode  Shape 
C=Circumferential  Mode  Shape 
D=Vertical  Responses 


16 


The  accelerometers  were  then  placed  at  the  coordinates 
required  for  the  first  test.  After  the  accelerometers  were 
installed  at  the  desired  location  and  orientation,  a  signal 
was  sent  to  the  accelerometers  to  insure  that  they  were 
functioning  properly,  and  that  all  cables  were  attached 
correctly.  The  low  pass  filter  was  then  set  at  the  desired 
frequency  (5,  10,  or  20  Hertz).  The  strip  chart  recorder  was 
utilized  to  display  the  accelerometer  output  while  the 
attenuation/gain  settings  on  the  SC-1  were  adjusted.  The 
gain  was  set  to  maximize  the  output  without  saturating  the 
magnetic  tape.  The  outputs  from  each  transducer  were  then 
recorded  for  up  to  twenty  minutes.  The  process  was  repeated 
for  low  pass  filter  settings  of  twenty,  ten,  and  five  Hertz. 

The  recording  of  signals  at  three  low  pass  filter 
settings  resulted  in  a  considerable  amount  of  data.  An  on 
site  spectrum  analysis  might  have  eliminated  two  of  the  three 
tests  since  preliminary  estimates  of  modal  frequencies  could 
be  obtained.  However,  a  spectrum  analyser  was  not  available 
for  use  on  site. 

Pre-printed  log  sheets  were  used  to  record  all  filter 
and  gain  settings  during  each  test.  These  were  referred  to 
during  the  tabulation  of  the  response  data  to  determine  the 
relative  magnitude  of  each  of  the  Fourier  spectra  generated 
by  the  spectrum  analyser. 

Figure  2.4  is  a  sample  of  the  pre-printed  log  sheets. 

Actual  log  sheets  may  be  found  in  Appendix  I. 
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TEST  LOG 


Date 
Time 


Test  No.  _ 
Rec.  Start 
Rec.  End 
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1 

2 

3 

4 
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Low  Pass 
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Attenuator 
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H-P  Filter 
I/O 

O 

as 

CO 

Sensitivity 
Volts 

Notes: 


Figure  2.4.  Pre-printed  Log  Sheet  Used  in  Ambient 
Vibration  Study 
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The  classification  of  each  test  is  also  noted  in  Table 
2.1.  Each  test  (excluding  the  calibration  test)  was  designed 
specifically  for  the  identification  of  the  predominant 
structural  frequencies  and  associated  mode  shapes. 
These  include  vertical  modes  as  well  as  radial  and  lateral 
modes . 

In  addition  to  the  18  tests  outlined  in  Table  2.1, 
free-field  tests  were  also  conducted.  This  was  accomplished 
by  placing  the  accelerometers  at  various  locations  on  the 
ground,  away  from  base  of  the  tower,  and  recording  the 
signals  in  the  manner  previously  described.  Free-field  data 
were  used  to  determine  if  any  predominant  site  frequencies 
exist  in  the  range  of  the  predominant  structural  response 
frequencies. 

Upon  completion  of  the  acquisition  and  recording  of  the 
data,  the  acceleration  records  were  analysed  at  Kansas  State 
University. 
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2.3  Data  Reduction 

The  spectrum  analyser  was  used  to  digitize  (sample)  the 
signals  and  to  convert  or  transform  them  to  the  frequency 
domain.  To  accomplish  this,  the  analyser  made  use  of  a  num- 
erical technique  called  the  Fast  Fourier  Transform  (22).  The 
Fourier  Transform  of  a  function  is  defined  as; 


oo 


F(t)  =  E   C  einut 


n 
n=-oo 


(1) 


where 


C  =  I 


r/ 


F(t)e-inUtdt 


(2) 


"The  Fast  Fourier  Transform  is  an  efficient  technique 

useful  for  a  computer  determination  of  the  response  in  the 

frequency  domain.  The  response  of  a  single  degree  of  freedom 

system  to  a  general  force  is  given  by; 

N-l  c  e2ir(nj/N) 

y(tj)  =  L       k(l-r*+2irn*  )  <3) 

n=0      n    n  n 

where  N  is  the  period,  the  damping  ratio  g  -  c//2km,  the 
frequency  ratio  rn  =U>n/cj  ,    and  the  coefficient 

N-l 

C-J  E  F(tj)e-2^^/N),  n=0>  1(  2 (N.n        (4) 

j=0 
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Either  equation  3  or  4  may  be  represented,  except  for 
sign,  by  the  exponential  function 

N-l 

A(j)=  £  A(0)(n)wJn,  (5) 

n=0 

where 

Wn  =  e2tfi/N  (6) 

The  evaluation  of  the  sum  will  be  most  efficient  if  the 
number  of  time  increments  N  into  which  the  period  T  is  div- 
ided is  a  power  of  2,  that  is, 

N  =  2M 
where  M  is  an  integer.  In  this  case,  the  integers  j  and  n 
can  be  expressed  in  binary  form.  For  the  purpose  of  illus- 
tration, consider  a  very  simple  case  where  the  load  period 
is  divided  into  only  eight  time  increments,  that  is,  N=8, 
M=3.  In  this  case,  the  indices  will  have  the  binary  repre- 
sentation     i^i^^+u^, 

n=n  +2n,+4n0, 
o   1   z 

and  equation  (5)  may  be  written  as 

1    1    1 
A(j)=  USE   A(0)(n)W8(J0+2Jl+4J2)(n0+2nl+4n2)    (7) 
n2=<^  ni"°  no"^ 

The  exponential  factor  may  be  written  as 

wJn=WS(Jln2+2J2n2+J2nl)W4n2J0W2nl(2Jl+J0)W^0(4J2+2Jl+J0)  (8) 
Note  that  the  first  factor  on  the  right  side  is  unity, 
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since  from  equation  6 

W8I  „e2iri(8/8)I  =CQS  2vl   +  1  sin2lrI  =  1( 

where  i=j-^n2  +  2j2n2  +  j 2ni  is  an  integer.  Therefore  only 

the  remaining  three  factors  need  be  considered  in  the  sum- 

ations .  These  summations  may  be  performed  conveniently  in 

sequence  by  introducing  a  new  notation  to  indicate  the 

successive  steps  in  the  summation  process.  Thus  the  first 

step  can  be  indicated  by 

1 

A(1)(j0,ni,n0)  =  £  A(0)(n2,n1,n0)W^n2J0) 
n2=0 

where  A   (no,n, ,Hq)  =  A    in  equation  5  .  Similarly,  the 
.  second  step  is 

1 
4(2)/,   ;  _  n  _  t-i  .(Dfj   _  „  u,2n-,  (21,+jn) 
A   (Jn,Jl'  o'    <L  A    (J()'  1'  Cr  8  ' 

n,=0 

and  the  third  step  (final  for  M=3)  is 

1 
A(3)(j0,j1,j2)  -S   A(2)(J0,J1,n0)wg0(4J2+2Jl+Jo). 
nQ=0 

( 3 ) 
The  final  result  A   (jn.ji.jo)  is  equal  to  A(j)  in  eq.5  . 

This  process,  indicated  for  N=8,  can'  readily  be  extended 

M 
for  any  integer  U=2  .  The  method  is  particularly  efficient 

because  the  results  of  one  step  are  immediately  used  in 

the  next  step,  thus  reducing  storage  requirements." 
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A  problem  associated  with  the  reconstruction  of  a  con- 
tinuous spectrum  by  discrete  data  is  known  as  aliasing.  This 
occurs  when  a  less-than-perfect  low  pass  filter  is  employed 
in  the  collection  of  data.  In  theory,  the  filter  should  re- 
move all  frequencies  above  the  limit  of  the  filter.  In  re- 
ality, however,  the  cut  off  rate  is  not  perfect.  This  allows 
the  reflection  of  response  peaks  above  the  cut  off  frequency 
to  be  displayed  below  the  cut  off  frequency.  As  an  example, 
consider  a  twenty  Hertz  sine  wave  sampled  at  a  rate  of  twenty 
samples  per  second.  There  will  be  one  sample  per  cycle,  the 
value  of  which  is  dependent  upon  the  time  the  sampling  be- 
gins. Because  the  sampling  rate  and  the  frequency  are  the 
same,  each  subsequent  sampling  will  have  a  value  equal  to 
the  first.  Figure  2.5  shows  a  sine  wave  with  a  sampling  rate 
equal  to  its  frequency.  Note  that  the  sampled  data  is  flat 
and  represents  a  zero  frequency  signal.  Thus  the  process  of 
sampling  data  at  the  same  rate  as  the  maximum  expected  fre- 
quency may  result  in  erroneous  reconstruction  of  the  data 
into  a  continuous  spectrum. 

Prior  to  the  analysis  of  the  signals,  several  combina- 
tions of  low  pass  filter  settings  and  sampling  rates  were 
tested,  and  it  was  determined  that  aliasing  would  not  be  a 
problem.  Therefore  the  signals  recorded  at  the  low  pass  fil- 
ter setting  of  twenty  Hertz  were  sampled  at  a  rate  of  40 
samples  per  second. 
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&   sampled  data 


Figure  2.5,  Twenty  Hertz  Sine  Wave  Sampled  at 
20  Samples  Per  Second 

A  total  record  of  1024  samples  was  used  per  recon- 
structed spectrum.  Thus,  25  seconds  of  response  data 
(1024/40)  were  used  to  calculate  each  spectrum.  A  total  of 
32  spectra  were  then  averaged  together  to  reduce  the  effect 
of  noise,  and  to  produce  the  Fourier  amplitude  plots.  Fig- 
ure 2.6  is  a  sample  of  the  spectral  amplitude,  as  a  function 
of  frequency,  for  the  signals  recorded  at  stations  5-0-N 
and  3-0-N. 

The  spectrum  analyser  was  also  used  to  determine  the 
relative  phase  of  any  two  responses,  as  a  function  of  fre- 
quency. In  addition,  the  coherence  function  was  also  deter- 
mined and  plotted.  The  coherence  function  (23)  is  given  as 

|Sxy(f) |2 


Kxy(f)   = 


Sx(f)Sy(f) 
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Figure  2.6.  Sample  Fourier  Spectra  From  Stations  5-0-N 
and  3-0-N 


where  x  and  y  represent  two  random  processes, 


Sxy(f)  = 
-co 

Sx(f)  = 


RxyWe"12*"  d-C, 


Rx(t)e-i21ff^d-Tr, 


Sy(f)  = 

-0O-' 


Ry(t)e"i2',ff'Cd'B'. 


The  autocorrelation  is  given  by; 

,T 


RxCC)  -  lim 
T-»oo 


TT~ 


x(t)x(t+tr)dt, 


Ry(t)  =  lim  ^jj- 


T-*CD 


y(t)y(t+T»)dt, 


and  the  cross-correlation  is  given  by; 

x(t)y(t+T»)dt. 


1 


RxyCE)  =  lim  j=r~ 
T-»eo 

-T- 


If  the  coherence  function  is  one  at  a  given  frequency, 
the  signals  are  highly  coherent,  and  probably  related.  If 
the  function  is  nearly  zero,  the  signals  are  highly  inco- 
herent, and  probably  unrelated. 
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Figure  2.7  is  a  sample  of  both  the  phase  and  coherence 
function  plots  from  the  reconstructed  signals  recorded  at 
stations  5-0-N  and  3-0-N. 

A  total  of  130  Fourier  amplitude  plots  and  65  relative 
phase  and  coherence  function  plots  were  generated  by  the 
spectrum  analyser.  These  may  be  found  in  Appendix  I. 
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Figure  2.7  Sample  Phase  and  Coherence  Function  Plots 
From  Stations  5-0-N  &  3-0-N 


III.  RESULTS 

Inspection  of  the  Fourier  amplitude  plots  revealed 
coherent  and  repeatable  peak  responses  at  seven  frequencies. 
These  frequencies  are  1.04,  1.44,  1.60,  1.84,  2.64,  2.88, 
and  3.52  Hertz.  The  spectral  amplitude  of  each  station,  at 
each  frequency,  was  measured  and  normalized  to  the  spectral 
amplitude  "of  the  reference  station  at  the  respective  freq- 
uency. These  normalized  values,  as  well  as  the  relative 
phase  values,  were  tabulated  and  may  be  found  in  Appendix  II. 

In  addition,  the  spectral  amplitude  of  the  reference 
station,  at  each  frequency,  was  normalized  to  the  spectral 
amplitude  of  the  reference  station  at  1.04  Hertz. 

The  tabulated  information  from  Appendix  II  was  also 
used  to  plot  lateral,  radial,  and  vertical  mode  shapes. 

Free  field  and  base  rocking  data  were  also  investigated 
and  modal  damping  ratios  were  calculated. 

3.1  Normalization  of  Reference  Spectral  Amplitudes 

The  spectral  amplitudes  of  the  reference  station,  at 
each  of  the  predominant  frequencies,  were  normalized  to  the 
amplitude  of  the  reference  station  at  1.04  Hertz.  This  pro- 
vided a  plotting  factor,  or  normalization  coefficient,  to  be 
applied  to  each  mode  shape  plot.  These  coefficients  may  be 
found  in  Table  3.1. 
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Frequency 

Measured  Response 

Normalization 

(Hertz) 

(MV) 

Coefficient 

1.04 

6.50 

1.00 

1.44 

3.90 

0.60 

1.60 

7.04 

1.08 

1.84 

29.68 

4.57 

2.64 

15.17 

2.33 

2.88 

8.23 

1.27 

3.52 

6.93 

1.07 

Table  3.1,  Reference  Responses  &  Normalization 
Coefficients 
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3.2  Mode  Shapes 

The  tabulated  response  and  relative  phase  information 
was  used  to  plot  lateral  (m)  ,  radial  (n)  ,  and  vertical  mode 
shapes,  at  each  of  the  predominant  frequencies. 

Lateral  mode  shapes  consist  of  responses  (normalized 
spectral  amplitudes)  normal  to  the  axis  of  revolution,  at 
a  specific  position,  and  at  varying  elevations.  Figure  3.2.1 
shows  typical  lateral  mode  shapes  of  a  cantilever  structure. 

Radial  mode  shapes  consist  of  responses  both  normal  to 
the  axis  of  revolution,  and  tangential  to  the  level  circular 
plane  of  the  structure,  at  varying  positions,  at  a  specific 
elevation.  It  must  be  noted,  however,  that  the  radial  mode 
shape  plots  are  defined  by  only  three  points,  and  therefore 
do  not  present  conclusive  evidence  concerning  the  radial 
mode  shape  value  of  n  at  any  level  or  frequency. 

Figure  3.2.2  shows  typical  radial  mode  shapes  for  an 
axisymetric  structure. 

Vertical  mode  shapes  were  plotted  and  examined  in  an 
effort  to  determine  the  contribution  of  base  rocking  to  the 
lateral  response  of  the  structure. 

When  more  than  one  value  is  available,  the  mode  shapes 
were  plotted  using  average  spectral  amplitudes  at  each 
station. 

The  lateral,  radial,  and  vertical  mode  shapes  will  be 
described  in  this  section. 
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n=0 


n=l 


n=2 


n=3 


Figure  3.2.2,  Typical  Circumferential  Mode  Shapes 
of  Axisymetric  Structures  (12) 
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1.04  Hertz 

The  lateral  mode  shape  of  the  tower  at  1.04  Hertz 
exhibits  a  maximum  response  at  level  3.  The  response  at 
level  3  is  1.50  times  greater  than  the  reference  response, 
and  is  out  of  phase. 

Levels  0  through  3.5  exhibit  responses  which  are  out 
of  phase  with  respect  to  the  reference  response,  while  the 
response  at  level  4  is  in  phase. 

Figure  3.2.3  is  the  lateral  mode  shape  at  1.04  Hertz. 
The  normalization  coefficient  is  1.00. 

The  radial  mode  shape  plots  at  both  levels  1  and  2  are 
indicative  of  an  n=0  mode,  however,  this  cannot  be  verified 
due  to  the  limited  number  of  data  points  available  for 
analysis. 

Figures  3.2.4,  and  3.2.5  are  the  radial  mode  shapes  at 
levels  1  and  2,  for  1.04  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.6. 
Levels  1  and  4  exhibit  upward  responses,  while  levels  0,2, 
and  3  exhibit  downward  responses.  The  maximum  vertical  re- 
sponse may  be  found  at  level  3,  and  is  0.37  times  greater 
than  the  reference  response. 
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1.44  Hertz 

The  lateral  mode  shape  of  the  tower  at  1.44  Hertz 
exhibits  a  maximum  response  at  level  4.  The  response  at 
level  4  is  1.87  times  greater  than  the  reference  response, 
and  is  in  phase. 

The  responses  at  all  levels  are  in  phase  with  the 
reference  response. 

Figure  3.2.7  is  a  plot  of  the  lateral  mode  shape  at 
1.44  Hertz.  The  normalization  coefficient  is  0.60. 

The  radial  mode  shape  plots  at  both  levels  1  and  2, 
are  indicative  of  an  n=2  mode,  however,  this  cannot 
be  verified  due  to  the  limited  number  of  data  points 
available  for  analysis. 

Figures  3.2.8  and  3.2.9  are  the  radial  mode  shapes  at 
levels  1  and  2,  at  1.44  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.10. 
Levels  0,  2,  3,  and  5  exhibit  upward  responses,  while 
level  1  exhibits  a  downward  response.  The  maximum  vertical 
response  is  found  at  the  top  of  the  tower. 
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1.60  Hertz 

The  lateral  mode  shape  of  the  tower  at  1.69  Hertz 
exhibits  a  maximum  response  at  level  4.  The  response  at 
level  4  is  1.50  times  greater  than  the  reference  response, 
and  is  in  phase. 

Levels  0  through  3  exhibit  responses  which  are  out  of 
phase  with  respect  to  the  reference  response,  while  the 
response  at  levels  3.5  and  4  are  in  phase. 

Figure  3.2.11  is  a  plot  of  the  lateral  mode  shape  of 
the  tower  at  1.60  Hertz.  The  normalization  coefficient 
is  1.08. 

The  radial  mode  shape  plot  at  level  1  is  indicative 
of  an  n=0  mode  while  the  plot  at  level  2  indicates  that 
n=l.  This  cannot  be  verified  due  to  the  limited  number 
of  data  points  available  for  analysis. 

Figures  3.2.12  and  3.2.13  are  the  radial  mode  shape 
plots  at  levels  1  and  2,  at  1.60  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.14. 
With  the  exception  of  level  5,  all  of  the  responses  are  in 
a  downward  direction,  and  from  levels  1  through  3,  seem  to 
increase  linearly.  The  maximum  vertical  response  occurs 
at  level  3. 
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1.84  Hertz 

The  lateral  mode  shape  of  the  tower  at  1.84  Hertz 
exhibits  a  maximum  response  at  level  5  (reference  station). 

Levels  3  and  3.5  exhibit  responses  which  are  out  of 
phase  with  respect  to  the  reference  response,  while  the 
responses  at  levels  1,  1.5,  2,  2.5,  and  4  are  in  phase. 

Figure  3.2.15  is  a  plot  of  the  lateral  mode  shape  of 
the  tower  at  1.84  Hertz.  The  normalization  coefficient 
is  4.57. 

The  radial  mode  shape  plots  at  both  levels  1  and  2 
are  indicative  of  an  n=0  mode,  however,  this  cannot  be 
verified  due  to  the  limited  number  of  data  points  available 
for  analysis. 

Figures  3.2.16  and  3.2.17  are  the  radial  mode  shapes 
at  levels  1  and  2,  at  1.84  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.18. 
With  the  exception  of  level  3,  all  of  the  responses  are 
in  an  upward  direction. 

The  maximum  vertical  response  may  be  found  to  occur  at 
the  top  of  the  tower. 
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2.64  Hertz 

The  lateral  mode  shape  of  the  tower  at  2.64  Hertz 
exhibits  a  maximum  response  at  level  5  (reference  station) . 

Levels  1.5,  3.5,  and  4  exhibit  responses  which  are 
out  of  phase  with  respect  to  the  reference  response,  while 
the  responses  at  levels  2,  2.5,  3  and  5  are  in  phase. 

Figure  3.2.19  is  the  lateral  mode  shape  at  2.64  Hertz. 
The  normalization  coefficient  is  2.33. 

The  radial  mode  shape  plots  at  both  levels  1  and  2 
are  indicative  of  an  n=0  mode,  however,  this  is  also  not 
verifiable. 

Figures  3.2.20  and  3.2.21  are  the  radial  mode  shapes 
at  2.64  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.22. 
Levels  2  through  5  exhibit  upward  responses,  while  level  1 
exhibits  a  downward  response. 

The  maximum  vertical  response  may  be  found  to  occur 
at  the  top  of  the  tower. 
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2.88  Hertz 

The  lateral  mode  shape  of  the  tower  at  2.88  Hertz 
exhibits  a  maximum  response  at  level  5  (reference  station) . 

Levels  1,  1.5,  2,  and  4  exhibit  responses  which  are 
out  of  phase  with  respect  to  the  reference  response,  while 
the  responses  at  levels  3  and  5  are  in  phase.  The  dashed 
line  on  the  lateral  mode  shape  plot  represents  the  uncer- 
tainty concerning  the  relative  phase  at  level  2.5. 

Figure  3.2.21  is  a  plot  of  the  lateral  mode  shape  of 
the  tower  at  2.88  Hertz.  The  normalization  coefficient  is 
1.27. 

The  radial  mode  shape  plot  at  level  1  is  indicative 
of  an  n=0  mode,  while  the  plot  at  level  2  seems  to  indicate 
that  n=3  or  more. 

Figures  3.2.22  and  3.2.23  are  the  radial  mode  shapes 
at  levels  1  and  2,  at  2.88  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.24. 
Levels  3  and  5  exhibit  downward  responses. 

The  maximum  vertical  response  may  be  found  to  occur 
at  the  top  of  the  tower. 
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3.52  Hertz 

The  lateral  mode  shape  of  the  tower  at  3.52  Hertz 
exhibits  a  maximum  response  at  level  5  (reference  station). 

Levels  2.5,  3,  and  4  exhibit  responses  which  are  out 
of  phase  with  respect  to  the  reference  response,  while  the 
response  of  levels  1,  1.5,  2,  3.5,  and  5  are  in  phase. 

Figure  3.2.27  is  a  plot  of  the  lateral  mode  shape 
of  the  tower  at  3.52  Hertz.  The  normalization  coefficient 
is  1.07. 

The  radial  mode  shape  of  the  tower  at  level  1  is  in- 
dicative of  an  n=l  mode,  while  the  radial  mode  shape  at 
level  2  seems  to  indicate  an  n=0  mode. 

Figures  3.2.28  and  3.2.29  are  the  radial  mode  shape 
plots  at  levels  1  and  2,  for  frequency  3.52  Hertz. 

The  vertical  mode  shape  may  be  found  in  Figure  3.2.30. 
Levels  1  and  5  exhibit  upward  responses,  while  levels  0, 
2,  and  3  exhibit  downward  responses. 

The  maximum  vertical  response  may  be  found  to  occur 
at  level  3. 


64 


CD 

CO 

CD 

3 

E 

*J 

4-J 

4J 

c 

G 

c 

3 

3 

G 

3 

c 

• 

M 

M 

M 

O 

O 

H 

O 

M 

t-H 

cu 

OS 

, 

oo 

*£> 

I— I 

r-* 

P-* 

M 

i — i 

O 

V 

m 

in 

in 

<f 

en 

CO 

r-. 

O 

> 

< 

o 

o 

o 

o 

o 

o 

O 

!-H 

en 

-3- 

o 

m 

09 

G^ 

CJ 

in 

TJ 

■ 

3 

0 

o 

•rl 

- 

fH 

CM 

r-» 

>— < 

m 

1 

<* 

*tf 

CO 

r-. 

o 

o 

o 

o 

« 

« 

„ 

„ 

CO 

l£ 

•4* 

r-. 

CM 

-—1 

r*. 

o 

CO 

\T\ 

^o 

■sf" 

<f 

CO 

o 

o 

o 

o 

o 

O 

o 

o 

o 

1—1 

H 

8J 

> 

J^ 

Jfi 

4P 

<U 

.— 1 

CM 

CM 

CO 

CO 

<r 

IT) 

J 

u 

M 

x 


a 
a. 

91 

■c 
5 
X 


W 
01 
U 
H 

J 


■H 


65 


o 

CO 


tu 

M 

ID 

J= 

cu 

4J 

4-1 

C 

3 

a 

a 

c 

• 

H 

Q 

5 

i — i 

H 

rH 

0) 

as 

01 

-a 

3 

0 

•H 

fH 

O 

r-* 

o 

iTt 

kO 

1 

lO 

m 

^D 

CM 

CM 

o 

o 

o 

O 

o 

s 

z 

0 

a 

1 

H 

•rl 

SB 

i 

o 

H 

4-) 

i 

o 

00 

O 

i 

4J 

C/3 

o 

ON 

r-l 

ON 

o 

^H 

* 

^H 

* 

^ 

2 


u 

a 


4J 

■e 


01 

■a 
o 
s 


■H 

■a 

5 

as 


CO 
CM 

CM 

CO 

91 

I 

■H 
[K 


66 


01 

a 

CS 

J3 

S 

4-J 

■U 

c 

0 

a 

9 

3 

• 

H 

M 

rH 

O 

O 

r-l 

a) 

os 

01 

■o 

3 

u 

■rl 

rH 

n 

■* 

i— 1 

*G 

in 

1 

to 

<r 

•* 

i— i 

CM 

o 

o 

o 

o 

o 

B 

35 

0 

2 

i 

H 

■rl 

SB 

I 

o 

H 

U 

l 

O 

00 

o 

a 

o 

o> 

, — 1 

c^ 

o 

4J 

t 

1 

1 

i 

w 

CNJ 

IN 

Csl 

CN 

eg 

67 


"Id 

U 

< 

< 

< 

< 

< 

tg 
V 

H 

a\ 

00 

vD 

<£ 

1 — 1 

a 

a 

a 

u 

5 

? 

.*. 

3 

•H 

a. 

0 

- 

o 

0 

a 

a 

a 

=: 

a 

a 

<ii 

T3 

3 

0 

CO 

t£ 

^H 

o 

o 

•H 

o 

.—1 

CM 

o 

CM 

i— t 

1 

o 

O 

O 

o 

o 

r-i 

OJ 

> 

Hi 

o 

r-H 

CM 

en 

LO 

o 
Q 


M 

CM 

in 


4-i 

1 


% 


u 


> 

o 


3.3  Free-Field  Tests 

Free-field  tests  were  conducted  at  various  locations 
away  from  the  base  of  the  tower.  These  locations  are  shown 
in  Figure  3.3. 

Figures  3.4,  3.5,  and  3.6  are  response  plots  generated 
from  the  three  free-field  tests.  Note  that  no  predominant 
frequencies  exist  in  the  range  of  the  predominant 
structural  response  frequencies.  This  precludes  the  poss- 
ibility that  the  predominant  structural  response  peaks  were 
artificially  induced  through  low  levels  of  seismic  activity. 

3.3  Base  Rocking 

Vertical  responses  at  the  base  of  the  structure  were 
examined  to  determine  if  any  base  rocking  occured  at  any  of 
the  predominant  frequencies. 

Table  3.4  summarizes  the  response  and  relative  phase 
information  from  test  17,  coordinates  0-0-V,  0-90-V, 
and  0-180-V.  Examination  of  Table  3.4  reveals  responses 
at  0-0-V  and  0-180-V  which  are  nearly  equal  in  magnitude 
at  all  frequencies.  With  the  exception  of  3.52  Hertz,  the 
relative  phase  is  also  the  same  at  each  frequency.  This  is 
not  indicative  of  a  significant  contribution  to  the  lateral 
response  by  base  rocking  about  the  90-270  axis. 
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Figure  3.3.   Approximate   Free-Field  Test   Locations 
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Station 

Frequency 

Response 

Phase 

o-o-v 

0.11 

In 

0-90-V 

1.04 

0.26 

In 

0-180-V 

0.11 

In 

0-0-V 

0.09 

In 

0-90-V 

1.44 

0.19 

In 

0-180-V 

0.09 

In 

0-0-V 

0.05 

In 

0-90-V 

1.60 

0.19 

In 

0-180-V 

0.05 

In 

0-0-V 

0.02 

In 

0-90-V 

1.84 

0.05 

In 

0-180-V 

0.02 

In 

0-0-V 

0.02 

In 

0-90-V 

2.64 

0.06 

Out 

0-180-V 

0.03 

In 

0-0-V 

0.05 

Out 

0-90-V 

2.88 

0.10 

Out 

0-180-V 

0.04 

Out 

0-0-V 

0.07 

In 

0-90-V 

3.52 

0.11 

Out 

0-180-V 

0.06 

Out 

Table  3.4.  Summary  of  Response  &  Relative 
Phase  Information  from  Test  17A 
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3.5  Modal  Damping  Ratios 

The  damping  ratio  may  be  computed  at  each  of  the 
predominant  frequencies,  provided  the  peaks  on  the  re- 
sponse spectra  may  be  isolated.  This  is  accomplished  by 
the  half-power  band  width  method.  The  method  states 
that: 

£= Ago/go  n 

where Ao5=  the  band  width  at  l/i/~2  of  the  peak  response 
at  an  isolated  frequency  (00n). 

Figure  3.7  shows  a  typical  response  plot,  complete 
with  half -power  measurements  and  calculations. 
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IV.  CONCLUSIONS 

An  ambient  vibration  study  was  conducted  on  the  hyper- 
bolic cooling  tower  at  Portland  General  Electric 's  Trojan 
Nuclear  Generating  Station,  near  Rainier,  Oregon. 

The  purpose  of  this  study  was  to  determine  the  dynamic 
characteristics  (natural  frequencies  and  associated  mode 
shapes)  of  the  tower.  The  mode  shapes  include  vertical  and 
radial,  as  well  as  lateral  responses. 

It  is  hoped  that  the  results  will  be  useful  in  the 
calibration  and/or  verification  of  results  obtained  from 
a  mathematical  model  (finite  element  analysis)  currently 
under  study  at  Kansas  State  University. 

Seven  predominant  frequencies  were  identified.  These 
are  1.04,  1.44,  1.60,  1.84,  2.64,  2.88,  and  3.52  Hertz. 
The  mode  shapes  associated  with  the  predominant  frequencies 
were  plotted. 

At  the  outset  of  this  research,  it  was  believed  that 
the  structure  would  be  found  to  behave  more  like  a  canti- 
levered  beam,  and  would  display  predominant  lateral  modes. 
Instead,  it  is  now  thought  that  the  predominant  modes  are 
circumferential  waves,  n=0,  1,  2,  3,  etc.,  each  with  a  set 
of  lateral  mode  shapes  unique  to  itself. 

In  as  much  as  there  are  insufficient  data  for  the  iso- 
lation of  n,  no  clear  conclusions  may  be  drawn  concerning 
the  relationship  of  m  to  n,  for  this  structure. 


77 


Examination  of  the  vertical  response  data  led  to  the 
conclusion  that  there  was  no  significant  contribution  to 
the  lateral  response  by  base  rocking. 

Free-field  data  indicated  that  no  predominant  site 
frequencies  exist  in  the  range  of  the  predominant  struc- 
tural response  frequencies.  This  precludes  the  possibility 
that  the  structural  response  peaks  were  artificially  in- 
duced by  low  level  seismic,  or  man  made  disturbances. 

The  principal  problem  encountered  with  the  testing  of 
a  structure  of  this  type  was  the  lack  of  accessibility  to 
the  various  areas  of  the  structure.  This  resulted  in  a 
lack  of  sufficient  data  for  the  identification  and  plotting 
of  radial  mode  shapes. 

Possible  remedies  include  the  construction  of  scaff- 
olding to  increase  the  accessibility,  the  use  of  optical 
methods  in  obtaining  data  (interferometry) ,  or  the  instru- 
menting of  a  tower  with  permanent  transducers,  during  the 
construction  of  the  tower.  This  would  supply  a  continuing 
record  for  analysis,  and  would  enable  the  researcher  to 
study  the  effects  of  aging  on  the  dynamic  characteristics 
of  the  structure. 
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Abstract 

A  series  of  ambient  vibration  tests  was  conducted  on 
the  hyperbolic  cooling  tower  at  the  Trojan  Nuclear  Gener- 
ating Station,  near  Rainier,  Oregon.   Acceleration  data 
were  collected,  and  analysed  via  the  Fast  Fourier  Trans- 
form.  Seven  predominant  frequencies  were  identified,  all 
below  five  Hertz,  and  associated  mode  shapes  were  plotted. 
These  include  vertical,  radial,  and  lateral  mode  shapes. 
Collection  of  data  for  identification  of  radial  modes 
was  hampered  by  inaccessibility  to  various  portions  of  the 
structure,  therefore,  accurate  representation  of  the  radial 
modes  was  not  possible.  Free-field  tests  were  completed. 
Modal  damping  ratios  were  calculated  for  three  frequencies. 
This  was  accomplished  using  the  Half-Power  Band  Width  meth- 
od. 


